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SUMMARY 


Approximate solutions are presented in this paper for three nonlinear ortho- 
tropic plate problems- The problems are to determine the deformations and forces 
in: (1) a thick plate attached to a pad having nonlinear material properties which, 

in turn, is attached to a substructure which is then deformed; (2) a long plate 
loaded in inplane longitudinal compression beyond its buckling load; and (3) a long 
plate loaded in inplane shear beyond its buckling load. For all three problems, 
the two-dimensional plate equations are reduced to one-dimensional equations in 
the y-direction by using a one-dimensional trigonometric approximation in the 
x-direction. Each problem uses different trigonometric terms. Solutions are 
obtained using an existing algorithm for simultaneous, first-order, nonlinear, ordi- 
nary differential equations subject to two-point boundary conditions- For each 
problem, reasons are discussed for choosing the trigonometric terms, ordinary differ- 
ential equations are derived to determine the variable coefficients of the trigono- 
metric terms, and sample results are presented. This solution procedure provides a 
quick, easy way to solve some nonlinear plate problems. 


INTRODUCTION 

Linear equations are often inadequate in the analysis of plates that appear in 
space vehicle structures and aircraft structures. Some problems of plate analysis 
involve stretching and bending coupling so that the strains depend nonlinearly on the 
deformations. Other problems may involve nonlinear material properties, such as 
those due to plasticity or nonhomogeneous materials- Strength or limits in deforma- 
tion may be considered in the design of such plates. Useful estimates of the defor- 
mations can usually be found from approximate solutions of the nonlinear equations of 
equilibrium. Meaningful strain distributions, however, require much more accurate 
solutions of the nonlinear equations. 

Many methods are available for solving nonlinear plate problems. A fairly com- 
plete survey, with examples of methods used before the high-speed computer was 
developed, is contained in reference 1. Analytical solutions of the differential 
equations are sometimes possible through special devices or inverse procedures. 

Direct procedures include iterative methods, such as Newton's method, and methods 
which make use of trial and error. For problems involving a parameter of small 
value, a regular perturbation method may be employed; for problems involving a bound- 
ary layer, methods of asymptotic integration may be employed. Weighted residual 
techniques are also available which include the energy method, the Galerkin method, 
collocation, and the method of least squares. Finite-element and finite-difference 
methods are often used now with the high-speed computer- These available methods 
have certain limitations, and solution to many nonlinear plate problems are not read- 
ily obtained by these methods. 

In the present paper, simplifying assumptions are used in a derivation from 
basic relations for plate problems to replace the nonlinear, partial (two- 
dimensional) differential equations of plate theory with nonlinear, ordinary (one- 
dimensional) differential equations. The derivation employs the principle of virtual 



work in conjunction with the assumption that the displacements may be represented by 
the first few terms of a Fourier series. Solution of the ordinary differential equa- 
tions, subject to the boundary conditions which arise naturally in the derivation, 
may be obtained by using the algorithm described in the next paragraph. 

An algorithm based on Newton* s method has been developed by Lentini and Pereyra 
(ref. 2) to solve a system of simultaneous first-order, nonlinear, ordinary differen- 
tial equations subject to two-point boundary conditions. The system of equations is 
of the form 


y’ = F(x,y) 


where y is the vector of dependent variables, x is the independent variable 
defined in the interval (a,b), and F is a vector, the components of which may be 
nonlinear functions of x and y. The boundary conditions of the problem are speci- 
fied by 

g[y(a),y(b)] = 0 

(The components of g identify each of the boundary conditions.) This algorithm 
uses finite differences with deferred corrections, and adaptive mesh spacings are 
automatically produced to detect and resolve mild boundary layers. 

Approximate solutions for three nonlinear plate problems are presented in this 
report. They are as follows: (1) a thick plate attached to a foundation having 

nonlinear material properties, which in turn is attached to a substructure, which is 
then deformed; (2) a long plate loaded in inplane compression beyond its buckling 
load; and (3) a long plate loaded in inplane shear beyond its buckling load. The 
first problem is important in designing thermal protection systems for space trans- 
portation vehicles and analyzing bonded joints. The second and third problems are 
important in structural design to take advantage of the stress-carrying ability of 
supported plates beyond the buckling load. The purposes of this paper are ( 1 ) to 
indicate the reasons for choosing the trigonometric Fourier series terms that are 
used, (2) to present the derivation of the differential equations that determine the 
coefficients of the Fourier series terms, and (3) to present sample results for each 
of the three problems. 


SYMBOLS 


44 55 66 


plate extensional stiffnesses 


A__ A A ^ A I- I- 

33p 44p 55p 


a,b 


pad extensional stiffnesses 
dimensions of rectangular plate parallel to x- and y-axes. 


C ,C ,C constant coefficients of cubic parts of stress-strain laws 
x y z 


respectively 
for pad 
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plate bending stiffnesses 


Dii,Di2/ 

^ 22'^66 


E,G 


E ,G 
P P 


M ,M ,M 
X y xy 

N ,N ,N 
K y xy 

N 

xav 

N 

xyav 

N 

zp 

T 

t 

u,v,w 

u,v,w 


Young’ s modulus and shear modulus for beam 
Young's modulus and shear modulus for pad 
length of beam 

bending moments in plate or beam 

inplate stress resultants in plate or beam 

average in plate 

average in plate 

normal stress resultant in pad 

transverse shear stress resultants in plate 

transverse shear stress resultants in pad 

thickness of pad 

thickness of plate or beam 

deformation of substructure in x-, y-, and z-directions 
displacements in x-, y-, and z-directions 


Ui,u2,Vq/Vi,v2, 

^x2 ' ^yO ' ^y1 ' ^y2 


arbitrary functions of y in equations (14), (15), and (30) for 
u, V, w, and 




“O'^OO 


cn 


u 


sh 


w 


x,y,z 


arbitrary functions of y in equations (41) for u, v, and w 

arbitary constants in ejqpressions (14) 
total axial displacement of beam 
applied compressive displacement 
applied shearing displacement 
total deflection normal to beam 
plate or beam coordinates 
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rotations in x- and y-directions 




^ > neutral surface strains 

^xy ' ^xz ' ^y z ) 


e ,y Y strains in pad 

zp xzp, yzp 

ic , K / K curvatures 

X y 

X half-wavelength in x-direction 

Prime (’) indicates differentiation with respect to y. Subscripts x and y which 
follow a comma indicate partial differential of the principal symbol with respect to 
the subscript. 


APPLICATIONS 

Derivations of ordinary differential equations are presented in this section for 
(1) a thick plate on a nonlinear foundation, (2) the postbuckling behavior of long 
plates in compression, and (3) the postbuckling behavior of long plates in shear. 

The choices of trigonometric Fourier series terms used for each plate problem, and 
the rationale for choosing these terms, are discussed. To clarify the derivation of 
the differential equations for the aforementioned problems, the governing equations 
for a simpler problem - the deep beam on a nonlinear foundation - is derived first. 
The deep-beam problem is the one-dimensional counterpart of the two-dimensional 
thick -plate problem and shows in a much simpler problem the steps involved in deriv- 
ing the differential equations of equilibrixam from the principle of virtual work. 


Deep Beam on Nonlinear Foundation 

As indicated in figure 1, the beam is attached to the top of a pad, which is 
attached to a substructure that may be deformed in the x- and z-directions (U and W). 
The effects of transverse shearing are included, in addition to extension and bend- 
ing. The pad resists transverse extension and shearing and has nonlinear material 
properties. The problem is to determine the strains in the beam and in the pad due 
to deformation of the substructure. Generally, the beam is considered to be much 
stiffer than the pad. 

Let u be the x-displacement of the neutral surface, w the deflection of the 
neutral surface, and the rotation. The neutral surface strain is then 


e = u* 

X 


( 1 ) 


the transverse shearing strain is 


‘xz 


= p + w' 


( 2 ) 
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and the curvature of the neutral surface is 


The bottom of the beam thickness t displaces u - (t/2)p^; thus, the rotation of 
the pad of thickness T is [u - (t/2)p^ - U]/T. The average deflection of the pad 
is (w + W)/2; the transverse shearing strain in the pad is then 


u - (t/2)p^ 
T 


U 


+ 


w' + W* 
2 


( 4 ) 


the sum of its rotation and slope. The beam is assumed to be much stiffer than the 
pad; accordingly, the slope of the beam w* is neglected in comparison with the 
substructure slope W* . Therefore, the transverse shearing strain in the pad is 


’xzp 


u - (t/2)p^ - U 

T ~2 


( 5 ) 


The extension of the pad w - W over the thickness T determines the transverse 
extensional strain in the pad as 


e 


zp 



( 6 ) 


The virtual work of the system for a beam of unit width and of length L is 


6n = r (N 6e + M 6 k + Q 6y + N 6e + Q 
XX X X X ' xz zp zp 


6y ) dx 
xp xzp 


( 7 ) 


where the coefficients of the virtual strains are force resultants. Substitution for 
the strains and integration by parts and rearrangement results in 


6n . JT t i e^) 6u + (a^ - m; - ^ a^) 

* {- 2 ; * f O «”] * =x ^"'0 

According to the principle of virtual work, the virtual and actual displacements 
must satisfy the same boundary conditions, but the virtual displacements are arbi- 
trary in the interior region. For equilibrium, the virtual work is an extremum; 


I 
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6n = 0. Therefore, equation (8) for virtual work can be replaced by the following 
differential equations and boundary condition terms: 


"N' + - Q =0 
X T xp 


Q Q =0 

^x X 2T 


-Q' + N =0 
X T zp 


(N 6u)^ =0 (M 

X 0 X '^x 0 


= 0 


«"'o = » 


( 9 ) 


Products of forces and virtual displacements appear in the boundary terms. The edges 
of the beam are free; therefore, = 0 at x = 0 and L. For other 

problems, the displacements can be taken to be equal to zero or combinations of 
forces and displacements (such that the above ejqpressions are satisfied) can be set 
equal to zero. 


A linear stress-strain law is assumed for the beam as follows: 


N = Ete = Etu' 
X X 


„ .si. 

X 12 X 12 ^x 


( 10 ) 




For the pad, a cubic nonlinear stress-strain law is assumed as follows (see fig. 1) : 
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N =ET(e -fCe^l 
zp P I zp z zp/ 

( Y + C ) 

^xzp x'xzp/ 


( 11 ) 


Q =GT|y +Cy 
xp p l^xzp x'xzp/ 


The equations expressing the stress-strain law assumed for the beam equations (10) 
and the equilibrium equations (9), together with the stress-strain law for the pad 
equations (5), (6), and (11) can be identified as six simultaneous first-order, 
ordinary, nonlinear differential equations with unknowns u, w, M^, 

and with boundary conditions. This system of equations may be readily solved by 

the algorithm discussed in the Introduction. The six equations in the desired form 
are as follows: 


u’ 



(12a) 


6 



(12b) 


12M 

3« =: 2 

Et 


W = 7;r - P 
Gt '^x 


N' = G 

X p 


u - (t/2)p - 


Gpt V" - (t/2)p - 

M' = Q ^ ^ 

X ’“x 2 


. E * C 

X p|_ T z\ T / J 


! . |l . .. [• ■ - ■ . |:]’l 

-•f- ■. [••'T-'" -f]’| 


(12c) 

(12d) 

(12e) 

(12f) 


Thick Plate on Nonlinear Foundation 

For a thick plate on a pad with nonlinear material properties and substructure, 
the idealization to a deep beam discussed in the preceding section would not be 
satisfactory if deformations across both the length and width of the plate are impor- 
tant. This derivation is the two-dimensional extension of the beam problem (see 
fig. 2 for the coordinate system) , and the substructure is permitted to deform V in 
the y-direction as well as the U and W in the x— and z— directions, respectively, 
for the beam. However, the U,V,W deformations considered are limited to 





ux 

2tdc 


u = Uqo + 

UqX + u 

i<y) 

sin — + V-y(y) 
a ^ 

sin 

a , 

I 

o 

> 

11 

> 



TOC 

2toc \ 


+ v^(y) 

COS 

— + ^^(y) cos 
a ^ 

a j 

) (13) 

W = Wg(y) 

+ W^(y) 


TlX 

2toc ] 

i 

cos 

— + Wo(y) cos 
a ^ 

a 



where the barred quantities are constants, and the indicated functions of y can be 
chosen as desired. As for the deep beam, transverse shearing is permitted in the 
thick plate, and the plate is considered to be much stiffer than the pad. The pad 
has nonlinear material properties and resists transverse extension and shearing in 
both the X- and y-directions. To make use of the algorithm described in the Intro- 
duction, it is necessary to make approximations and derive ordinary differential 
equations to replace the partial differential equations of plate theory. 


In some linear problems solved by the use of Fourier series, the loading terms 
are expanded in Fourier series, and the solution is in the form of a Fourier series 
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which is exact term by term. For various approximations/ the series may be limited 
to the first few terms of the series. A similar approximate method is used here for 
a nonlinear problem. The substructure deformations represented by equations (13) can 
be considered to be loading terms / and/ after identifying the first few terms of the 
series for the solution/ no further terms are considered (and no further approxima- 
tions are made ) . 

Trigonometric terms of the same form as those considered for U/ V, and W in 
equations (13) are assumed for the plate displacements n, V/ and w in the x-/ y-/ 
and z-directions as follows: 





TlX 


2tcx 

u = Uqo + 

UqX + U. 

i(y) 

sin — + u^(y) 
a ^ 

sin 

a 




7DC 

2tix 


O 

> 

II 

> 

+ v-|(y) 

cos 

— + v^(y) cos 
a ^ 

a 





TlX 

2m 


II 

o 

+ w^(y) 

cos 

— + wo(y) cos 
a ^ 

a 



(14) 


Similar trigonometric terms are assumed for the plate rotations |3^ and p as 
follows: 


nx 2 idc 

Px = sin — + sin 


7TX 

Py = + Pyl(y) COS — + 


27DC 

a 


(15) 


The neutral surface strains and curvatures in the plate are 


^X ^/X 


Yxy = ^^,y + ^x 


fC-- p 


X/X 




= V 




*^XZ Px ^,x 


<y Py,y 


(16) 


Yyz = py + -,y 


'‘xy 


Px,y Py,x 


The commma before a subscript denotes partial differential with respect to that sub- 
script. The strains in the pad are given as 


u - (t/2)pj^ - U W 

Y = + 

’xzp T 


/X 

2 


(17a) 


8 



2 


(17b) 


V - (t/2)p - V W 

Y = ^ + 

yzp T 


e 

zp 


w - W 


T 


By siibstitution from equations (14) and (15), the strains and the curvatures are 


n TUK 2% 2tix 

— + u^ — cos 

^ ^ 'a a ^ a a 


■Jtx 2 ti:x 

e = Vq + vij cos — ^2 cos 

•* a a 



^yz 


7TX 


2ii;x 


PyO + + (Pyl + '^i) (Py2 '"z) ~ 


K 


X 


Px1 


71 TlX 

- COS — + 

a a 


Px2 


271 2tix 

— cos 

a a 


•^y = PyO + Py1 P + Py2 ^ 

/ 71 \ 7TX / 2ti \ 27tX 

<xy = (Pil - Py1 r (Px2 - Py2 r ) T 


(17c) 


(18a) 


(18b) 


(18c) 


(18d) 


(18e) 


( 18 f ) 


(18g) 


(18h) 
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'xzp 






2 Pxl - 



sin 


7tx 

a 


+ 




^x2 


T 


2 a y 


sin 


27TX 

a 



e 


zp 



w. 


- W, 


cos 


nx 



cos 


27tx 

a 


With no loss in generality, the constants appearing in u can be set equal to 
corresponding constants appearing in U as follows: 


*^00 “ ^00 


'^0 “ ^0 


The stress-strain laws for the plate are 


^x “ ^11^x ^12^y 


\ " ^22^ ^12^x 


2x ^44^xz 
2y “ ^55^yz 


Mx D-| + ^12*^y 


\ “ ^22*^y ^12*^x 


^xy ^GG^xy 


^xy ^GG^xy 


The stress-strain laws for the pad are 


2xp ^44p^^xzp ^x^xzp^ 


2yp ^SSp^Yyzp “*■ ^yYyzp^ 


^zp ^33p^^zp **“ ^z^zp^ 


(18i) 


(18j) 


(18k) 


the 


(19) 


( 20 ) 
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The form of the force and moment resultants for the plate in terms of the trigonomet- 
ric terms in the x-direction with coefficients/ functions of y, is similar to the 
form of the strain for the plate: 


Ttx 2 toc 

a a 


TDC 2tdc 

Ny = NyQ(y) + Ny^(y) cos — + Ny 2 (y) cos 


TUX 2ii;x 

^xy = N^yl(y) sin — + ~ 


= 


M,, = 


ux 2'fvx 

Mxo(y) + Mjjl(y) cos — + M 2<y) oos 

a a 


TlX 2 tix 

Mvn(y) + M„i(y) cos — + M 2(y) oos 

a -* a 


Y “ ^ yi 


IDC 2 TOC 

^*xy = Mxyl(y> sin — + M^y 2 <y) sin — 


Qx = Qxi < y ) r ^ 


TtX 2 TCX 

Qy = Qyo(y) + Qyl(y) cos — + Qy2<y) cos 


However, the form of the force resultants for the pad include terms 


TOC 


2tix 


2 ^ = 2^1 <y> sin — + Q-(y) sin + ... + Q„„^(y) sin 


. Gtix 




“"xp ^xpl 
3 =0 

yp ypo 

^ = N ^ 

zp zpO 


a ^xp2 


xp6 


2_ = 2_„(y) + 2_,(y) cos + Qyp 2 <y> ^ + 2.._.(y) cos 


"ypi 

zp1 


TOC 2 tdc 

*=°s — + N^p2^^^ cos — + ... + N_^(y) cos 


®yp6' 

*zp6< 


6toc 

a 

6tix 




( 21 ) 


( 22 ) 
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The first few of these coefficients are 


2xp1 = ^44 p[yxzp1 - I ylzp^ YxzplYxzp2^j 

2xp2 “ ^44p^Yxzp2 “ 2 *^x ^Y^zplYxzp2 2’ Yxzp2^j 

[ /a 3 2 1 2.32 \1 

YyzpO ^ylYyzpO 2 YyzpOYyzpl 2 YyzpOYyzp2 4 YyzplYyzp2 1 1 


^SSpjjyzpl ^^y ^YyzpOYyzpl YyzpoYyzplYyzp2 4 Yyzp 

^55p^^yzp2 ^^y^YyzpoYyzp2 2 YyzpOYyzpl ■*■ 4 Yyzp2^^ 


1 3 . J. 2 ^ 

■4 Yyzp1 ^ YyzplYyzp2 


[ / 3 3 2 3 2^32 V] 

^zO ^z^^zO 2 ^z0^z1 2 ^z0^z2 4 ^z1^z2jj 

^33p|^^z1 ^ ^^z^^zO^zl ^ ^z0^z1^z2 ^ 4 ^z1 * 4 ^z1^z2^J 

^33p|^^z2 ^^z^^z0^z2 2 ^z0^z1 4 


where 


,1 - (“1 - 1 - ” 1 )/'' - r s 


All the other pad strain coefficients are similar. (See eqs. (18).) 
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The virtual work of the system for a plate of width b and length a is 


6n 



( N 6e + N 8e + N 
XX y y xy 


6v + Q 6v 

xy X xz 





+ M 

6k + M 

6k 

+ M 

6k 

X 

X y 

y 

xy 

xy 

+ Q 

6y + 

Q 

6y 

+ N 

xp 

xzp 

yp 

yzp 

zp 



dx dy 


Substituting for the force and moment resultants from equations (21) and equa- 
tions (22), substituting for the strains from equations (18), and integrating 
X results in 



+ 


+ 


+ 




J['' [\l I ♦ “X2 «"2 ¥■ * ♦ "yl + •■ya 

Vl (®"i ■ ‘’t i) * '*>W2(^"2 - ^’2 ?) 2*1 

Q o(6p ^ - 6w^ — I + 2Q ^(6p „ + 6w') + Q , + 6w* ) 

x2 V ^x2 2 ^ ® yi yi 1 

0 y 2 '« Py 2 ^ ®" 2 > + “*1 « P *1 I * "*2 ‘^*2 

^y1 ^^y1 * ”y2 *^y2 ”*yl(^^x1 " ^^y1 “xy2^^^x2 

Py2 ¥} * e*pl f - I «P*l) * =xp2 ?(‘»2 - I ^^* 2 ) 



* ^2ypo t(^ 0 “ 2 *^yoj 2yp1 “ 2 *^ylj 

* «yp2 fK - I ^^, 2 ) ^ I ‘“0 ^ Vl f 

* ^2 5 **2] 


( 24 ) 


over 


(25) 
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Because of the integration properties of products of orthogonal functions, only 
the first few coefficients for the forces in the pad appear. Integration by parts 
results in 


6n 




, + N , -21 

xy2 x2 a 


6u + (“2N* + — + /-N' . - N . — + 6v, 

2\v0 T/0 


y1 xy1 a T 
2N 


+ /-N* -N — + 6v + /-2Q' + — Sw 

I y2 xy2 a T/2ly0 TJO 

/ ^zpl^ ^ 2% ^zp2\ 

I ^y1 ^x1 a T / 1 y ^y2 x2 a ^ / 

^ (-%1 “xl i * 2x1 - %’-) «Pxl 


6w, 


- + M — 
xy2 x2 a 


+ Q, 


, 66 + f-m' +20 y£2_\ 58 + 

x 2 2 T / '^x 2 I yO '^yO T / '^yO 


2xp2^ 


+ Q. 


- 66 + f-M' - M + Q 

y1 2T / Py1 ^ y2 xy2 a ^y2 2T 


Q -.t 

yp2 


/-M*. - 

I 

) J 


M ^ - 
xy 1 a 

dy 


-6u.+N . 6u + 2N 6v^ + N ^ 6v- + N^ 6v^ + 2Q ^ 6w^ 

2 \ xyl 1 xy2 2 yO 0 y1 1 y2 2 ^yO 0 

+ Q ^ Sw + Q ^ 6w^ + M ftp + M ^ 6P ^ + 2M ^ ftp ^ + M ftp ^ 
y1 1 y2 2 xy1 x1 xy2 x2 yO yO y1 ^y1 

+ M ftp 

Y2 ^y2yo 


Thus, the principle of virtual work requires satisfaction of the following differen- 
tial equations and choice of boundary conditions: 


xy 1 




N 


xy2 





^ypO 

T 


(26a) 


(26b) 


(26c) 
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(26d) 


N* = -li N . + 


■^YP1 


y1 a xyl T 


N* = - 15 N + 

“y2 a ”xy2 


(26e) 


= I ”xl ^ 2x1 - 


2xpl^ 

2T 


(26f) 


M* = — M + Q - ■ 

xy2 a x2 ^2 2T 


(26g) 


M' n = Q. 


2yp0^ 


yO *y0 2T 


(26h) 


M* = - - M . + 
y1 a 


xyl 2yi 


2yp1^ 

2T 


(261) 


M* = - — M 

+ 0 - ^yp^^ 

(26j) 

y2 a xy2 

«y2 2T 

Q. = 

^yO T 


(26k) 

Q* 1 - -Q . ^ 

y1 x1 a T 

(261) 

Q'o = -Q o — + 
y2 ^x2 a 

^zp2 

T 

(26m) 

N , = 0 

xy1 110 

N ^ 6u„|^ = 0 
xy2 2 1 0 

(26n) 

N 6v = 0 

yO 0 lo 

N . 6v = 0 N ^ 6v P = 0 

y1 1 |0 y2 2 10 

(26o) 

M , 6p = 0 

xyl ^x1 10 

M o Sp ^ = 0 

xy2 ^x2 |0 

(26p) 
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6 p = 0 

yO ! 0 

M , 6 p = 0 

yi yl |0 

M ^ 66 ^ 
y 2 '^y 2 

O 

tl 

Xi o 

( 26 q) 

O 

11 

o 

Q 1 

yi 1 


Q „ 6 w 
y 2 2 

0 = ° 

( 26 r) 


The boundary condition assumed for the results presented in this paper is that 
the edges are free. Therefore/ the stress and moment resultants are zero; so at 
y = 0 and b 


N = N ^ = N ^ 
xyl xy2 yO 


N 


yi 


= N 


y2 


= Q..n = Q... = Q. 


"yo 


yi 


■^y2 


M = M ^ 
xyl xy2 


0 


M 


yO 


= M 


yi 


= M 


y2 


0 


The system of first-order, ordinary differential equations to be solved for this 
problem are presented in terms of the following 26 unknowns: 


r ^^2' "o' 


^1' 


V 


/W # w^/ w / 

2' o' 1 2 


Pxl' 


x2 


^yO' ^yl' 


y2 


xyl' \y2' "yO' "y1 ' V' V' ''yi' '^y2' ^'xyl' ‘'xy2' ‘'yO' "yi' "y2 


N 


N 


N 


M 


M 


M 


M 


M 


Equations (26), which were obtained from the virtual work, present 13 of the differ- 
ential equations used. The remaining 13 equations were obtained from the strain- 
displacement relations (eqs. (18)), the stress-strain relations (eqs. (19)) for the 
plate, and equations (21). These remaining equations are as follows: 


u 




(27a) 


u: 



N 


xy2 


66 


(27b) 


V' 

o 




(27c) 
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(27d) 


V 


1 




u 


1 


V- = ^ - (hl\^ U 


(27e) 


®yO 

w* = - 6 

0 A '^■ 

55 


yO 


(27f) 


w 


= fi 

1 - Py1 


(27g) 


w 


B 

2 A53 Py2 


(27h) 


M 


S' = + JE R 

^x1 Dgg a Py1 


(27i) 


M 


6' = ^ fi 

Px2 -^a Py2 


(27j) 


yO D__ 
22 


(27k) 


”v1 

S' . = 

^y1 D 

22 



8 

,D22^ ^X1 


(271) 


**y2 

P' « = - 

y2 

22 



^ li2E s 

(022 ^ Px2 


(27m) 


17 



In addition, the stress-strain laws for the plate directly result in 
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and, after some elementary manipulation, result in 
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Sample results for this problem are presented subsequently. 
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Long Plate in Compression Loaded Beyond Buckling 

The analysis of the postbuckling behavior of long, rectangular orthotropic 
plates in longitudinal compression is presented in this section. The plate has a 
length a and a width b, and the long edges are supported. The following sketch 
of a buckled plate identifies the periodic deflection w of half -wavelength X: 
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At buckling, the equations for w are satisfied by 


/ , . *nx 

w^(y) sxn — 


and the equations for u and v are satisfied by ^ ^o(b " 

respectively, where u^^ and Vq are constants. The study of the postbuckling 
behavior of a plate presented in reference 3, which solves for the unknowns in the 
order of their importance (perturbation method) , indicates that the next group of 
terms to be considered are 


The next 


u^Cy) sin 

V 2 (y) cos 
term to be 
w^(y) sin 


2itx 

\ 


2tpc 

X 

considered is 

3tcx 

X 


However, for up to three or four times the buckling load, the coefficient w^ was 
found to be very small compared with w^ for a long plate (ref. 3). Accordingly, 
the assumed displacements are 




2-rtx 

T' 


2-rex 


V = Vp(y) + v^(y) cos ^ 


TCX 

W = W^ (y) sin — 


(30) 


Thus, the deflection w is assumed to be sinusoidally periodic with ha If -wave length 
The deflection w is exact at buckling, the displacements u and v are sinu- 
soidally periodic with ha If -wave length X/2, and u has an extra term which is 
linear in the x-direction associated with the constant u^^^ which is specified. No 
further approximations are made. 


A Levy -type solution is developed for the postbuckling behavior of long plates 
in compression. The neutral surface strains and curvatures as given by Von Karman 
nonlinear plate theory are 
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From the stress-strain law for an orthotropic plate, the form of the stress resul- 
tants can be identified as follows: 
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where 
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The virtual work of the system is 
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Substitution of equations (31) and (32) into equation (34) and integration over x 
results in 
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Integration of equation (35) by parts results in 
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where/ by definition. 
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Thus, the principle of virtual work requires satisfaction of the following 
ential equations and choice of boundary conditions: 
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Since the variations must satisfy the same boundary conditions as the actual dis- 
placements, the variation signs were omitted in the boundary terms in equations (39). 
For the y boundaries, the simply supported, straightedge boundary condition 
requires that at y = 0 and b 


w =N =u = v = M =0 
1 yO 2 2 y1 


The clamped, straightedge boundary condition requires that at y = 0 and b 


w 


= N 


yO 


= V = P = 0 


Thus, NyQ = 0 throughout the region for both boundary conditions considered. The 
simultaneous first-order differential equations to be solved can be written in terrnis 
of the unknowns 
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with Vq determined by integration after the other unknowns are determined. The 
following differential equations for these unknowns have been determined by manipula- 
tion of equations already stated: 


N 


= 


xy2 

^66 


2n 

X 


1 71 o 

^2 


N o 
22 










1 


w; = 


P 


N 


xy2 


= 2u|il2 

^ h 2 



(40a) 


(40b) 


(40c) 


(40d) 


(40e) 


23 


(40f) 
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Equations (40) are the eight ordinary differential equations to be solved for this 
problem- Sample results for this case are presented subsequently. 


Long Plate in Shear Beyond Buckling Load 

Analysis of the postbuckling behavior of long, rectangular orthotropic plates in 
inplane shear is presented in this section- The plate has a width b, and the long 
edges are supported. The following sketch of a buckled plate identifies the skewed 
periodic deflection w of ha If -wave length X: 



“sh/2 


At buckling (see ref. 4), the equations for w are satisfied by 

TlX TOC 

w (y) sin t — + w (y) cos t — 
s X c ^ X 

The subscripts s and c refer to sine and cosine. The equations for u and v 

are satisfied by u , ) and 0, respectively. The choice of trigonometric terms 

shyb 2/ 

for the shear problem is not as straightforward as for the compression problem. 

Guided by experimental and theoretical results for shear which indicate that the 
buckle mode does not change much in the initial postbuckling range (as with 
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compression) , the postbuckling mode in shear is given by the buckling mode and a few 
elementary terms beyond. Accordingly, the assumed displacements for this problem are 


. . ... 2-iix , ^ 2 tix 


V = v-(y) + V (y) sin + v (y) cos 

US A. C A. 
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(41) 


No further assumptions are made. Following steps similar to those for compression 
loading/ the equations of equilibrium as obtained from virtual work are as follows: 
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Differential equations obtained from the stress-strain law are as follows; 
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By definition. 
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Additional equations from the stress-strain law are as follows: 



(45) 
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Equations (42), (43), and (44) (except for the fourth equation of (42) and the fourth 
equation of (43)) are the 18 ordinary differential equations to be solved for this 
problem. The fourth equation of (42) shows that N^q is a constant or zero and the 
fourth equation of (43) permits a solution for Vq when N^q has been determined 
and when the rest of the equations have been solved. The applied shearing displace- 
ment Ugj^ enters the problem through boundary condition on Uq« Sample results for 
this case are presented subsequently. 

SAMPLE RESULTS AND DISCUSSION 

Sample results are presented for each of the three plate problems. The solution 
to each of the problems appears to be more efficient by this method than any other 
method, including finite elements and finite differences. 
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Thick Plate on Nonlinear Foundation 


The thick plate considered has dimensions of 5 in. by 5 in. by 2 in., with 
extensional and flexural properties as follows: 


Ai 1 = A ^2 = 51 674 Ib/in. 

^12 ” Ib/in. 

= A^^ = 21 000 Ib/in. 
44 55 66 


'11 


12 


66 


^22 ” in-lb 

3100 in-lb 
7000 in-lb 


The pad is 0.16 in. thick with properties in the initial (linear) range given 
by; transverse extension modulus E = 500 psi and transverse shear modulus 
G-| =G 2 = 12.5 psi. Therefore, 


33p 


44p 

The substructure 


80 Ib/in. 
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55p 

is deformed such that 
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0.0003 


(i--) 
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0.12 sin ^ + 0.03 cos ^ cos — + 0.01 
b b a 
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sxn -T — cos 

b a 


where U, V, and W are measured in inches. 

The solution was obtained by the present method for the thick plate on the pad 
and deformed substructure. The values of the moment M and stress resultant 
at the center of the plate for various values of the pad cubic stiffness coefficients 
are shown in figure 3. When the coefficients are zero, the solution of the equations 
is exact. The moment doubles in the range shown as the cubic stiffness coefficients 
increase. The stress resultant increases by a factor of about 9 for the same range. 


Long Plate in Compression Loaded Beyond Buckling 

Equations of equilibrium have been derived in the section entitled "Applica- 
tions” which may be solved directly by using the algorithm described in refer- 
ence 2. For a given value of the applied compressive displacement u^n' 
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prescribed values of the dimensions, the material properties, and the ha If -wave length 
X, the equations may be solved and the average load intensity ^xav 
determined as follows: 
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xav 
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2b X 



2X 


N 

X 


dx dy 


For finite plates, equilibrixjm paths can be determined in the postbuckling range 
which may be associated with buckling modes or changes in mode. In addition, from 
examination of the energy, these paths can be identified as being stable or unstable. 
For an infinitely long plate, consideration may be restricted to wavelength, and, by 
comparison with results for finite plates, the wavelength of the stable path corre- 
sponds to the wavelength that gives minimum energy. An example of a path which is 
not of interest beyond buckling is one that gives zero deflections. The load of 
interest for postbuckling studies is the one that corresponds to an applied displace- 
ment that is in the postbuckling range and on the equilibrium path corresponding to 
the wavelength of interest. For this problem, these conditions are satisfied if 
(1) the applied displacement is larger than its buckling value, (2) the results give 
nonzero deflections, and (3) various half-wavelengths are tried until a minimum 
energy is obtained. 

Characteristic curves are plotted in figure 4 for the compression of long, iso- 
tropic plates and long ±45® laminated composite plates with quasi-isotropic lay up. 
The average stress intensity coefficient is plotted as a function of the applied 
displacement coefficient for the long edges held straight (inplane) and simply 
supported or clamped. The isotropic curves apply to isotropic metal or composites 
with an isotropic lay up. The ±45® laminate curves apply to graphite -epoxy filamen- 
tary material with properties given by the following dimensionless quantities: 
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-0.431 


For the isotropic plate, both these quantities are unity, and for both the isotropic 
and ±45® laminate. 


A 
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22^11 

D 

11 22 
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The slope of the load-displacement curve is a measure of the overall plate 
stiffness. As shown in figure 4, this curve is a straight line with a slope equal 
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to one prior to buckling. After buckling/ this line changes slope according to the 

Di2 + ^22^11 

boundary conditions and according to — — and r • (See ref. 5.) The 

Y^ir 22 ^11^22 

±45® laminate has a slightly lower postbuckling slope than the isotropic and, there- 
fore, has a slightly lower postbuckling stiffness. 


The only assumptions made in this analysis are that Von Karman theory applies 
and that the deflection w and the displacements u and v are sinusoidally 
periodic with half -wavelengths \ and X/2, respectively. Comparisons with experi- 
ment show that Von Karman theory is satisfactory. The deflection pattern is exact at 
buckling, and, by comparison with other work, this method should be accurate to three 
times the buckling load or higher. The results given herein for simply supported 
edges are the same as those presented in reference 6. The results given herein for 
clamped edges have not been presented previously. The curves for the clamped edges 
are similar to the curves for the simply supported edges, except for higher slopes in 
the postbuckling range. This difference indicates that clamped plates are stiff er in 
the postbuckling range. 


Long Plate in Shear Loaded Beyond Buckling 

For the plate loaded in shear, equations of equilibriiim have been derived in the 
section entitled ''Applications*' which may be solved directly by using the algorithm 
described in reference 2. For a given value of the applied shearing displacement 
Ush' prescribed values of the dimensions, the material properties, and the 

half-wavelength X, the equations may be solved and the average shear load intensity 
may be determined as follows: 

1 

N = II N dx dy 

xyav 2b\ xy 


As with the long plate loaded in compression, to get a postbuckling solution, (1) the 
applied displacement must be larger than its buckling value, (2) the results must be 
on the equilibrium path where the deflections are nonzero, and (3) various wave- 
lengths must be tried until the wavelength is obtained which corresponds to minimum 
energy. Characteristic curves are plotted in figure 5 for the shearing of long iso- 
tropic plates and long ±45® laminated composite plates with quasi-isotropic lay up. 
The average stress intensity coefficient is plotted as a function of the applied 
displacement coefficient for the long edges held straight (inplane) and simply sup- 
ported, or clamped. In addition to the parameters mentioned for compression (see 

^ 11^22 " ^12 " ^^ 12^66 

ref. 6), the shear case depends on the value of rr -- . Comparing 
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1^22 


the postbuckling slopes shown in figures 4 and 5 shows that the stiffness of plates 
buckled in shear is higher than the stiffness of plates buckled in compression. 


Again, the only assumptions made in this analysis are that Von Karman theory 
applies and that the deflection w and the displacements u and v are sinusoi- 
dally periodic with half-wavelengths X and X/2, respectively. Comparisons with 
e>q)eriment show that Von Karman theory is satisfactory for shear loading. The 
deflection pattern and the deformation pattern are skewed, unlike the compression 
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case, and the deflection pattern is exact at buckling. The method should also be 
accurate to at least three times the buckling load. The results given herein for 
simply supported edges are the same as those presented in reference 6. The results 
given herein for clamped edges have not been presented previously. As with the 
plates loaded in compression, the curves for clamped edges are similar to the curves 
for simply supported edges, except for higher slopes in the postbuckling range. The 
higher slopes indicate that clamped plates are stiff er in the postbuckling range. 


CONCLUDING REMARKS 

A solution procedure is presented in this paper for three nonlinear plate prob- 
lems. This procedure is demonstrated by deriving one-dimensional equations in the 
y-direction using a one-dimensional trigonometric approximation in the x-directipn 
for three different (two-dimensional) plate problems and by presenting sample results 
for these problems. For a plate on a nonlinear foundation, the solution is a non- 
linear extension of a Fourier series solution of a linear differential equation in 
which the loading terms are also expanded in a Fourier series to get an exact solu- 
tion, term by term. The postbuckling problems are solved by a nonlinear extension of 
the (linear) Levy “type solution for long plates. The trigonometric terms are exact 
for the linear range of stiffness for the plate on the foundation. The trigonometric 
terms are exact for initial buckling in the postbuckling problems. The sample- 
problem results for the plate on a nonlinear foundation show that with an increase in 
nonlinear stiffness of the foundation, the neutral-surface direct stress resultant is 
affected much more than the corresponding moment. Postbuckling plate problems 
studied include isotropic plates and a ±45® filamentary composite laminated plate. 

The ±45® laminated plate is not as stiff as the isotropic plate in the postbuckling 
range for the compression or shear loadings. There is not much difference for com- 
pression, but they differ considerably for shear. For both plates and both loadings, 
the postbuckling stiffness is higher for clamped edges than it is for simply sup- 
ported edges. 
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Figure 1.- Analytical model of deep beam attached to a substructure by a pad with 

nonlinear properties. 
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Figure 4*- Characteristic curves for plates in compression. 
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Figure 5.- Characteristic curves for plates in shear* 
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